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SUP CONDITIONS WITH WALL CATALYS IS AND RADIATION 
FOR MULT I COMPONENT, NONEQUI LI BRIUM GAS FLOW 

I. INTRODUCTION 


To accurately predict the heat transfer to the Space Shuttle during high 
altitude flight, the proper slip conditions at the gas -wall interface must be 
employed. In those cases in which the mean free path is small in comparison 
with the characteristic length of the body (Kn « 1) , the boundary condition at 
the surface of a solid body is that the gas velocity relative to the body is zero 
and that the temperatures of the body and gas are identical. However, these 
conditions are approximate and, when the Knudsen number is finite, the bound- 
ary conditions must be refined. Specifically, the density at the wall is suffi- 
ciently rarefied that the velocity, temperature, pressure, and concentration at 
the wall are different from that of the immediately adjacent gas. These "slip" 
or "jump" conditions will be derived. 

At high altitude-low Reynolds number flight, the Knudsen number is 
finite and the continuum model of the gas breaks down in regions of large 
gradients such as near a cold body. Hence, the Navier-Stokes description 
cannot be used in the region of gas a distance on the order of one mean free 
path away from the wall ( 1] . This region of gas is commonly referred to as 
the Knudsen layer. 

Since there are not sufficient collisions between gas molecules for the 
Navier-Stokes approximation to be valid in the Knudsen layer, the more 
rigorous Boltzmann equation should be utilized. The boundary conditions for the 
Boltzmann equation would be such that the solution would match the solution of 
the Navier-Stokes equations in the bulk flow and satisfy the kinetic boundary 
conditions at the wall. 

Without actually solving the complex Boltzmann equation in the Knudsen 
layer, Patterson [2] and Shidlovskiy 1 3] have used a slip model to find an 
approximate Navier-Stokes solution to the flow in the knudsen layer. This 
approach will be used herein. Scott 1 4] extended the approach of Shidlovskiy 
to include a multicomponent mixture with diffusion and wall -catalyzed atom 



recombination; however, second order terms in the velocity and temperature 
slip equations were neglected and it contained errors. The complete expressions 
with the effects of radiation are included herein, and simplified results for 
engineering applications are discussed. 

Sample calculations are made with single and multicomponent species 
gas slip conditions to show the necessity of including the multicomponent 
species in the slip conditions for catalytic and noncatalytic surfaces. 


II. ANALYSIS 


By matching the species, momenta, and energy fluxes at the outer edge 
of the Knudsen layer to the difference between the incident and reflected fluxes 
at the wall, the jump in properties at the wall are obtained. The fluxes are 
assumed to be constant across the Knudsen layer. These fluxes are calculated 
from moments of the distribution function. The distribution function comes 
from a Chapman-Enskog expansion for a multicomponent mixture obtained 
through a variational method. The Chapman-Enskog expansion, which for 
small departures from local translational equilibrium leads from the more 
exact Boltzmann equation to the approximate Nuvier-Stokes equation, is used 
to give solutions for the slip conditions compatible to the order of appi'oximation 
of the Navier-Stokes equations . 

It is assumed that the molecules specularly reflected from the wall are 
associated with the slip temperature while the molecules diffusely reflected 
from the wall are associated with the wall temperature. That is, the specularly 
reflected molecules have no accommodation with the wall while the diffusely 
reflected molecules are fully accommodated. The accommodation coefficient 
for specie "i", 0., is the fraction of those gas molecules (atoms) colliding with 

the surface that stick and subsequently leave the surface with the energy asso- 
ciated with the wall temperature. Some of the atoms which stick recombine 
with other atoms to form molecules. The recombination coefficient for specie 
"i", y . , is the fraction of those incident atoms that recombine at the surface. 

Changes in the internal energy of reflecting molecules are neglected. 

Little is known about the actual complex mechanism of gas interaction 
with the surface of a solid body. Accurate measurements of accommodation 
coefficients are not available because the physical process is not the same for 
different materials. For instance, surface adsorption of the gas can be a 
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factor or the molecules may condense on the surface and then evaporate . The 
physical shape of a surface can also affect the gas-wall interaction. When the 
surface is smooth the molecules tend to reflect specularly; however, in reality, 
the surface is more likely to be rough, having interstices in which the molecules 
may be temporarily trapped. Reflection of the molecules from a surface of 
this type is considered diffuse. The molecules strike the boundary with com- 
plete loss of the tangential mass velocity and escape after achieving partial or 
complete accommodation with the wall, the orientation of emission being ran- 
dom. The molecules might also react with other atoms or molecules on the 
surface before being emitted. Thus until more specific knowledge of gas -wall 
interaction is obtained, the aforementioned model will be utilized. 

Using the notation shown in Figure 1 and denoting 0(v) as a convected 
property such as mass, momentum, or energy, the incident flux of that prop- 
erty normal to the wall is 


V > 0 00 

F * = f f f V 0(vj f (V) d 3 V 

J J J y g 

-60 -60 -00 J 


(1) 


where f (V) is the velocity distribution function at the edge of the Knudsen 
s 

layer . The specularly reflected flux is proportional to 


F» 


oo oo 0O 

f f f v 0(v) f (V , -V , V ) d 3 V 
j j j y sx y z 
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The diffusely reflected flux is proportional to 


F = 
w 


f 

-60 



( 3 ) 


where 


f 

w 


(V) is the velocity distribution function at the wall. 




Figure i. The Knudsen layer showing general fluxes and coordinate 
axes. The temperature as a function of normal distance 
is schematically overlayed. 

.The general flux balance for the convected property is 

Z F i = Z F i ,+ Z (1 ■ V F i ,+ 

i i i 


ZK 


V 


w. 

1 


( 4 ) 


where the sum over species is not required for the species mass balance. The 
first term on the right is the incident flux, the second term is the specularly 
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reflecting flux (incident minus those that stick) , and the third term is the 
diffusely reflecting flux (those that accommodate to the wall minus those that 
recombine) . 

Equation (4) must be altered if radiation is included. If radiation 
stresses and radiation energy density are neglected but one-dimensional radia- 
tion heat flux is included, then only the energy flux balance differs from the 
form of equation (4) . Since the gas is rarefied, an optically thin or emission 
dominated gas is assumed. The Knudsen layer is assumed transparent because 
of its thickness. This approximation is consistent with assuming constant 
fluxes across the Knudsen layer. The energy flux balance with radiation is then 


Ye. = V E.' + E * + Y (1 - 0.) E.» + (1 - e) E » + Y (0. -y.)F. W 

U l u i R v li R V l l l 


+ E. 
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where 


E * = - E J - f 4o <7 T 4 dy 
R R •' n J 


(5a) 


y=6 


and 


W A 

E - - e ff T 4 
R w 


(5b) 


The velocity distribution function for a multicomponent mixture per- 
turbed out of equilibrium is 


f. (V) = f °(V) [1 + ip.(V)] 

l i l 


( 6 ) 


where 



Vin T - B :W n + 

lO 0 



d. 
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(Oa) 
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and f °(v) is the Maxwelliam velocity distribution for the ith species. The 
coefficients in equation (6a) are functions of the dimensionless velocity, 


W. 

1 



(6b) 


( i ) 

and contain the constants a . „ , a , , b.„, c., which are constants determined 

iO it iO iO 

from the Curtiss and Hirschfelder I 5] variation problem in the first approxima- 
tion for a mixture. These constants are related to the transport properties and 
functions of the collision integrals 1 6] . The vector d. is related to the diffu- 
sion velocity of the jth species and is defined as •* 



£_ 

m. 

J 


X 


- V n. X. 
j V i ij 


(6c) 


where X. are external forces. The total mass-average velocity, which is to 
be used later, is 


= 1 V n. 
P V ] 


m v 

j j 


(7) 


where v. is the sum of the bulk velocity v^ and the thermal velocity V. for 
species j . 

By inserting equations (1), (2) , and (3) into equations (4) and (5) with 
the use of equations (6) and carrying out the triple integrations, the equations 
for the slip conditions are obtained. All accommodation coefficients 0 are 

assumed equal to 0 . The slip properties are related to gradients at the outer 
edge of the Knudsen layer, the properties at the wall and, if radiation is 
included, the temperature profile. 
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Concentration slip: 



Pressure slip: 



( 9 ) 

Velocity slip: 
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(concluded) 


III. BOUNDARY CONDITIONS 


By equating the net mass flux of atoms to the surface to the rate of con- 
sumption of atoms at the wall from surface recombination, Scott [4] established 
the wall recombination rate constant in terms of the recombination coefficient 
and the thermal velocity as follows: 



Thus the thermal velocity at the wall limits the maximum rate constant for a 
fully catalytic (y. - 1) wall. However the rate constant for a fully catalytic 

wall is often assumed, for simplicity, to approach infinity. 

The net mass flux-rate of comsumption balance just described can be 
solved to obtain 



( 14a) 


9 



where V . is the diffusion velocity defined as 


CO 30 OO 

n S V . = f f f V . 4>(V) f. (V) d 3 V 

1 VI J J J yi IS 

v -OC -OC -00 


( 14b) 


Equations (8) through (12) and (14) form a coupled set of slip conditions 
to be utilized with flow field calculations. They are applicable to a multicompo- 
nent, nonequilibrium gas reacting catalytically on a surface with radiation and 
diffusion. 


IV. SIMPLIFICATION 


The constants a. y , a^, la ^ and c. y J in the slip conditions, equations 

(8) through (12) and (14) , can be determined from the collision integrals by a 
variational technique in which they are solutions to sets of simultaneous equa- 
tions [5] . Due to the complexity of the collision integrals, a great simplifica- 
tion results if one relates these constants directly to the transport properties 
[5]. This simplification amounts to using: (1) the relation for a i() 



n m 

_i l 
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(15a) 


T 

where D is the thermal diffusion coefficient, (2) the relation for a.. 

i il 



(15b) 


where K is the thermal conductivity, (o) the relation for b 



i 





(15c) 
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where ») is the viscosity coefficient, and (4) the relation for c 


D.. 

i] 


lz£L c (j) 

2 nm. v m. iO 
J i 


( 15d) 


where D.. is the diffusion coefficient. The diffusion velocity can also be 
ij 

written in terms of transport properties as 


V 


y i 



2 


m. 


D.. d. 

i] jy 


1 d T DinT 

n. m. i dv 
ii J 


(16) 


Substituting equations (15) and (16) into the slip conditions yields the slip con- 
ditions in terms of transport properties. Second order terms in the velocity 
and temperature slip expressions are now dropped. 

Concentration slip: 



Pressure slip: 
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I 



12 
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Temperature slip: 



Before any further simplifications are made, it is interesting to note 
what effects influence the magnitude of the slip conditions. The concentration 
slip, from equation (22) , is directly proportional to the net normal mass flux 
of species i and inversely proportional to the thermal slip velocity of species 
i . The velocity slip is related to the tangential shear stress, plus a thermal 
creep term over species i due to thermal diffusion and heat conduction in the 
tangential direction, minus a diffusion term in the tangential direction, all 
divided by the pressure over the thermal slip velocity. The pressure slip is 
proportional to the pressure at the wall, minus the normal shear stress, plus 
a term containing normal thermal diffusion and normal thermal conductivity, 
all divided by one plus a normal diffusion term . The temperature slip is pro- 
portional to the diffusely reflected energy flux, minus the normal flux of 
energy (including radiation) , all divided by a term proportional to the density 
times the thermal slip velocity. 
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Further simplifications can be made by making the assumptions consis- 
tent with Fick's law of diffusion: (i) thermal diffusion is negligible, (2) diffu- 
sion due to pressure gradients and external forces is neglected, and (3) a 

binary mixture of atoms and molecules is assumed so that D.. = D,„. By 

ij 12 

assuming that the normal shear stress for species i is the same as that of the 
bulk mixture (r . = r or n.r). = mi) , the concentration slip and temperature 

yyi yy i 1 

slip are simplified. Bj' assuming that the energy conducted for species i is 
the same as that of the bulk mixture, the velocity slip and pressure slip ai’e 
simplified. The above assumptions do not significantly influence the slip 
conditions . 

With these assumptions equations (17) through (22) become 
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and 



These expressions simplify to those presented by Shidlovskiy [ 3] for a 
single species if one assumes that the jump conditions are small and radiation 
is neglected. These jump conditions with radiation for a single species gas are 
as follows: 
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where an elastic sphere molecular model has been assumed such that 


ip /i — A- 
8 V2 ^/RT 


and 
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15 k 
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(34 a) 


(34b) 


Note from equation (33) that the radiation being absorbed at the wall causes 

T to be larger than that without gas emission. This term can also be inter- 
s 

preted in effect as an increase in the wall temperature. 


V. SLIP CONDITIONS NEAR A STAGNATION POINT 

Using spherical coordinates, nondimensionalizing by defining u - u/v^ , 

v = v/v^ , p - P/P^ , T = T/T 0 , 7} - v/Vo , P = P/pJJ and r = ~ r/r h ’ and 

rearranging terms, the expressions needed to determine slip conditions near a 
stagnation point for a chemically reacting gas are as follows: 
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(39) 
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Equations (35) through (39) with the equation of state form a complete 
set of boundary conditions for a multicomponent, nonequilibrium gas reacting 
catalytically on a surface. Knowing the wall temperature and thus the recom- 
bination rate at the wall, the recombination coefficient (y. ) can be obtained 

from equation (13) . The accomodation coefficient, (), can be assumed equal to 
one until better experimental values are obtained. 

For a completely non catalytic wall, y equals zero and equation (36) 

reduces to the radial gradient of the mass fraction of species "i" equal to zero. 
For a completely catalytic wall, y^ equals one. Using values of y^ of one and 

zero gives simplified expressions which show the effect of catalyticity . 


VI. SAMPLE CALCULATIONS 


To determine whether or not there is any significant difference in using 
the single species gas slip conditions or the multicomponent reacting gas slip 
conditions presented here, an accurate stagnation point solution is utilized. The 
full Navier-Stokes equations are solved from the free stream through the merged 
shock layer to the body using the concept of local similarity [7] . A rarefied 
portion of the Space Shuttle orbiter reentry trajectory with catalytic and non- 
catalytic walls is used as an example. The flight conditions are shown in 
Table 1 . 

Since the pressure and mass fraction at the wall are not specified A 
PRIORI, equations (35) and (37) will not be utilized. This is equivalent to 
assuming that the pressure slip is the same as the pressure at the wall and that 
the mass fraction slip is the same as the mass fraction at the wall. 

Table 2 presents a summary of the effect of using the more realistic 
multicomponent gas slip conditions. Results are presented for a spherical nose 
the same size as the Space Shuttle orbiter and for a typical reentry flight 
condition. 
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It is evident from Table 2 that erroneous results would be obtained using 

the single species slip expressions in lieu of the multicomponent expressions . 

Using the multicomponent species boundary conditions, the velocity slip is only 

52 percent as much and the temperature slip is only 33 percent as much as that 

using the single species boundary conditions at the highest altitude considered 

for either type wall. The heat transfer coefficient rises the most (14 percent) 

at 96. 2 km using the multicomponent species expressions and a noncatalytic 

wall. The heat transfer coefficient continues to rise with altitude over the 

single species expressions with the use of the more exact boundary conditions 

for a catalytic wall. For the catalytic and noncatalytic wall conditions, the 

multicomponent species slip conditions show a reduction of approximately 15 

percent in the skin friction coefficient at the highest altitude. From Table 2 it 

is evident that the use of the more exact multicomponent slip conditions causes 

a rise of approximately 15 percent in the heat transfer coefficient and a drop of 

approximately 15 percent in the skin friction coefficient at the highest altitudes 

when compared with results using the single species slip conditions. Other 

trajectories may well show even stronger variations in C.. and C . 

hi 1' 

Since this is not necessarily the most severe example , it appears that 
the more complex multicomponent reacting gas slip conditions must be used for 
accurate prediction of heat transfer and skin friction coefficients. 
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TABLE 1. FLOW FIELD FREE STREAM PROPERTIES 
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TABLE 2. EFFECTS OF SINGLE SPECIES AND MULTICOMPONENT SPECIES 
SLIP EQUATIONS FOR ORBITER TRAJECTORY 


o 

Single 

Species 

d o 3. O 

3 3 3 3 
05 CD CD 05 
t- CO CO O 
’T CO 

O h •- N 

n c c 

too 
3 3 3 'X 
05 *05 CO -3 
— LO — 05 
3 CO lO CO 
1- C- 3 rr 
CO 3 O Cl 

o 

3 m 

1 O CD 

-30 
— 0 3. 
o CO 

C: o 

3 3 3 3 

05 05 K 05 
O OO [' C5 
L3 O L5 X 
LO Tj 1 N O 
M C ® 

c t o 

CD CO CD >3 
Cl GO 1- 05 

DO CO rr Cl 
C ’f L' X 
CO 3 75 3 

o“ 

tn 
o o 

it « 
3 O 
.3 Q. 
CO CO 

O CO O 13 

co *r o - 
r- i' co a 

—< CO LO 3 

n CO Ol CO 

rT 3 co 

I - Cl 1- — 
Cl T LO I s " 

O 

S 05 
1 O O 

5 £o 

-GO 

— o a 

O CO 

t % C5 Ci 
Cl CO O I- 
r- o o i- 

•TH TT w |- 

efj 3 3 -V 

■h C IC5 « 
1 - rr Cl OQ 
ci rr 3 i" 

To 

H 

05 

q> O 

*3) 

.s g. 

CO CO 

C5 rr Cl — « 
IO 1- O f 
CDClOl- 

— * c> o* tj» 

iO CO Cl -4 
— LO 
3 oo 35 3 

*< Cl W 'l’ 

5 

3 05 

1 c o 

5 o 

3 c o 
— o a 

CO O CO 

ci i- o oi 

— <M 

i- o no in 
-I i>i ^ 

i - co io co 

Cl CO CV 1- 
3 33 1 - to 

^ ^ ^ ^ 

K 

05 
<U CJ 

*3) o 

3 £ 
CO CO 

3 3 

3 3 3 <£, 
co CD .3 .3 

O 1- CD 05 
r /> iO IM rr 
O O CO — < 
PO I* H TJ- 

O O -rH 

0 3. 

1 1 I = 

3 1 - CD W 
Tf 3 | - 'O 
O c h c 
CO 3 -h rr 
O O ~ 

"co 

3 

3 

O 

3 05 

1 o O 

i £ o 

— o a. 
^0(0 

'i. <t> 'x. 

3 3 3c: 
CD CD CD *05 
O h Cl 1- 
CO 1- O rr, 
Cl C H P5 
CM iO l- i- 

o o o o 

<X/ 'X. 

3 3 3 3 

CD CD CD CD 
Cl 3 C5 1- 
O 1- CO 35 
o o -« CO 
CO If; 1- 1 ~ 

o o o o 


I . Noncatalytlc Wall 
Altitude (km) 

87.43 

91.50 

90.2 

100. 

11. Catalytic Wall 

Altitude (km) 

87 . 43 
91.50 
80.2 
1 00 . 


23 




REFERENCES 


1. Kogan, M. N. , Rarefied Gas Dynamics , Plenum Press, New York, 1969, 
pp. 367-400. 

2. Patterson, G. N. , Molecular Flow of Gases , John Wiley and Sons, New 
York, 1956, pp. 119-125. 

3. Shidlovskiy, V. P., Introduction to the Dynamics of Rarefied Gases , 
American Elsevier Publishing Co. , New York, 1967, pp. 60-68. 

4. Scott, C. D., "Viscous Reacting Flows with Wall Slip and Catalysis 
Applied to Spheres in Arc Jets and Flight," AIAA Paper No. 74-174, 

1974. 

5. Hirschfelder, J. O., Curtiss, C. F. , and Bird, R. B., Molecular 
Theory of Gases and Liquids , John Wiley and Sons, New York, 1964, 
pp. 468-490. 

6. Scott, C. D. , "Wall Boundary Equations with Slip and Catalysis for 
Multicomponent Nonequilibrium Gas Flow," NASA TM X-58111, 1974. 

7. Kumar, A. , "Hypersonic Rarefied Transition on Blunt and Slender 
Bodies, " Ph.D. dissertation at Indian Institute of Technology Kanpur, 
1973. 


24 



I 


APPROVAL 


SLIP CONDITIONS WITH WALL CATALYS IS AND RADIATION 
FOR MULTICOMPONENT, NONEQUILIBRIUM GAS FLOW 

By Dr. William L. Hendricks 

The information in this report has been reviewed for security classifi- 
cation. Review of any information concerning Department of Defense or Atomic 
Energy Commission programs has been made by the MSFC Security Classifica- 
tion Officer. This report, in its entirety, has been detei’mined to be unclassi- 
fied. 

This document has also been reviewed and approved for technical 
accuracy. 



Homer B. Wilson, Jr. 

Chief, Thermal Environment Branch 



Werner K„ Dahm 


Chief, Aerophysics Division 

S~\ // ' if 

( ) /'I farm 0 * -*' 

r . i* ■ - t - — 

J.'-'A,, Lovingood 


Director, Systems Dynamics Laboratory 


25 



